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Abstract 


Cumulus clouds are observed to produce a very small fractional cloudiness, 
yet in the GLAS climate model they are assumed to fill an entire grid area 
(~500 km X 400 km). To study the effect of such excessive cumulus cloudiness 
on the radiational heating and, consequently, on other aspects of the model 
climatology, we have performed an experiment, for a July simulation, in which 
the cumulus cloudiness is neglected completely for purposes of the solar and 
terrestrial radiation parameterizations. The results are compared with those 
of a control run, in which 100% cumulus cloud cover is assumed. 

The net solar radiation input into the earth-atmosphere system is more 
realistic in the experiment, and the model's underprediction (by about 40 W m“2) 
of the global-mean outgoing thermal radiation at the top of the atmosphere is 
somewhat reduced. 

Comparison of the experiment and control suggests that there exists a 
positive feedback between cumulus convection and the radiation field. Treating 
the cumulus clouds as having zero cloudiness produces fewer deep cumulus clouds; 
apparently, the radiative effects of cumulus clouds tend to favor further con- 
vective activity. 

The upper troposphere is warmer in the experiment, especially in the middle 
latitudes of the Northern (summer) Hemisphere. The vertical shear of the zonal 
wind is correspondingly reduced. The surface air temperature increases over 
land, and the thermal lows over the continents intensify. 

The precipitation maximum over the tropical oceans shifts southward in the 


experiment. 


1 . Introduction 


A number of cloud-radiation studies have previously been made with the GLAS 
(Goddard Laboratory for Atmospheric Sciences) climate model (Harman et al . , 1980; 
Shukla and Sud, 1981); these studies ware reviewed by Herman (1980). In a "trans- 
parent cloud" experiment, Herman et al. (1980) eliminated the cloud interactions 
first with solar radiation and then with terrestrial radiation. They concluded 
that the treatment of all clouds as transparent to the radiation stream increases 
the global-mean radiative flux into the top of the atmosphere by 16 W m~2, in a 
"fixed cloud" experiment (Shukla and Sud, 1981), the spatial distribution of 
clouds was computed once and fixed in time. The results showed that a fixed 
cloud distribution acts as a zonally asymmetric, stationary thermal forcing, 
and produces a significant change in the large-scale dynamical circulation. 

These studies demonstrated the importance of cloud-radiation interactions 
for climate modeling, but they did not distinguish between cumulus cloudiness 
and supersaturation cloudiness. This distinction is important, first because 
of the difference in the fractional cloudiness associated with these two cloud 
types, second because of the different geographical and seasonal distributions 
of cumulus and stratus clouds, and finally because of the different vertical 
extents of these two cloud forms. 

In the GLAS climate model, simulated clouds of any kind are assumed to fill 
an entire grid area, which is typically about 500 km x 400 km in the horizontal 
and 1 - 2 km in the vertical. The simulated cloud cover is used to determine 
the radiative heating and cooling of the atmosphere-earth system. In nature, 
the fractional cloudiness produced by cumulus convection is observed to be 
small (Malkus, 1958; Krishnamurti, 1968). This suggests the need for a parame- 
terization of subgrid-scale convective cloudiness for use in the radiation 
parameterizations of the climate model. 


The theory of subgrid-scale fractional cloudiness is still in an early 
stage of development. The interaction of finite clouds with solar radiation 
has been studied by McKee and Cox (1974, 1976), Davies (1978), and Davis et 
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al . (1979a,b), and the interaction of finite clouds with infrared radiation has 
been discussed by Liou and Ou (1979), and Harshvardhan and Weinman (1982). A 
theory to determine the subgrid-scale cloudiness associated with shallow cumuli 
was proposed by Albrecht et al. (1979) and Albrecht (1979, 1981), while subgrid 
scale stratiform cloudiness has been parameterized by Sundqvist (1978, 1981). 
These studies represent progress towards the parameterization of cloudiness on 
unresolved scales, but the problem is so vast that much additional work is 
needed before a general parameterization suitable for global climate models can 
be created. In particular, there currently exists no theory to determine the 
fractional cloudiness associated with deep cumuli and their associated anvils 
and cirrus blowoff. Much of the cloudiness in the ITCZ and over the tropical 
and summer-hemisphere continents is believed to be associated with such convec- 
tive complexes (e.g. , Houze, 1977; Ogura and Chen, 1977). 

Lacking a physically based parameterization to determine the subgrid-scale 
fractional cloudiness associated with cumulus convection, it may be better to 
neglect the cumulus cloudiness altogether than to treat it as overcast. The 
purpose of the present study is to examine the impact on the simulated general 
circulation as the cumulus cloud cover is changed from 100% (the control) to 
zero (the experiment) for purposes of calculating the solar and terrestrial 
radiative transfer. As will be shown, this change leads to some improvement in 
the simulated radiation budget of the earth. It also reveals the strength and 
nature of the sensitivity of the model to fractional cumulus cloudiness. It 
thus provides guidance and motivation for future efforts to incorporate into the 
model a realistic parameterization of subgrid-scale cumulus cloudiness. 
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2.^ Model description and initial conditions 

The GLAS climate model has evolved from the GISS GCM, described by Somer- 
ville et al . (1974) and Stone et al « (1977), which in turn, was based on the 
3-level 1970 version of the GCM developed by A. Arakawa and Y, Mintz at UCLA. 

The modal uses a form of the a coordinate system (Phillips, 1957), with nine 
layers, all of the same a thickness, and a 4 x 5 degree latitude -longitude 
grid, which is modified in the polar regions as discussed below. The variables 
are staggered in the horizontal according to the B-grid scheme of Arakawa and 
Lamb (1977), and in the vertical according to the scheme of Lorenz (1960). The 
upper boundary of the model is at 10 mb. The prognostic atmospheric variables 
are the surface pressure, the zonal and meridional components of the horizontal 
wind, the temperature, and the water vapor mixing ratio. The prognostic boundary 
variables are the bulk ground temperature and ground wetness, and the snow depth. 

Although the current model is similar to the GISS GCM presented by Somer- 
ville et al. (1974) and Stone et al . (1977), significant changes have been made 
in both the finite-difference schemes and the physical parameterizations. As 
described by Halem et al. (1979) and Herman and Johnson (1978), the model incor- 
porates a "split grid," in which the number of grid points on a latitude circle 
is systematically reduced near the poles. This allows the use of a ten-minute 
time step, with only weak longitudinal smoothing at high latitudes. The idea 
is similar in principle to that proposed by Kurihara (1965), but is applied at 
only a few interfaces separating latitudinal bands within which the horizontal 
resolution is uniform. The resolution in each latitudinal band is twice that 
of its poleward adjacent neighbor. The differencing scheme for horizontal 
advection maintains the quasi-conservative properties of kinetic energy and mean- 
square vorticity, as in the Arakawa scheme. Although written to handle an 
arbitrary number of bands, the geometry used in the present model has five 
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latitudinal bands, consisting of a central region extending from -66" to +66", 
bounded by one set of bands f^oro 70" to 78®, and a second set from 82" to 86"i 
The latitudinal resolution is 4" in all bands. The longitudinal resolutions 
are 5®, 10", and 20", respectively. 

Once every simulated half hour, a sixteenth-order Shapiro filter (Shapiro, 
1970) is applied, in the longitudinal direction, to the sea level pressure, and 
to the potential temperature and the wind components on the o surfaces. The 
filter was introduced as an ad hoc device to suppress the "checkerboard" noise 
in the prognostic fields, which arises from the inability of the B-grid to sim- 
ulate the geos trophic adjustment process at the two grid interval level (Arakawa 
and Lamb, 1977). 

The planetary boundary layer (PBL) parameterization was formulated by 
Katayama and modified by Somerville, et al. (1974) and Sud and Abeles (1981). 

The model incorporates a prognostic soil moisture, with parameterizations 
of runoff, snow melt, and potential evapotranspiration, as developed by Lin et 
al. (1978). 

The Matsuno forward-backward time differencing scheme is used. This scheme 
tends to damp high frequencies, and so it helps to control computational noise 
in the model, but at the cost of almost doubling the computing time required 
for the dynamics. 

We use the cumulus parameterization developed by Arakawa (1969) for the 
three-level UCLA GCM; but as modified for use in a nine-level GCM by Somerville 
et al. (1974; see also Helfand, 1979). Although this Arakawa cumulus parameter- 
ization predates the Arakav;a-Schubert cumulus parameterization (Arakawa and 
Schubert, 1974; Lord and Arakawa, 1 981 ; Lord, 1982; Lord, Chao, and Arakawa, 
1982), many key concepts employed by Arakawa and Schubert were already present 
in the 1969 paper, including the cumulus mass flux, a spectrum of clouds, and 
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closure formulated in terms of the stabilization of the environment by the 
clouds. Helfand (1979) modified the parameterization to include very deep 
clouds and cumulus friction, but these modifications are nob incorporated into 
the present model. 

The model also includes latent heat release due to large-scale saturation, 
which occurs when the relative hv’*midity exceeds 100%. 

The short-wave radiation parameterization is that presented by Laois and 
Hansen (1974). It includes absorption by ozone, water vapor, and clouds. 

The long-wave radiation parameterization of the current model is based on 
the method of Wu (1976), Wu et al. (1978), and Wu (1980), as described by 
Krishnamurthy (1982). The parameterization includes a water vapor transmittance 
that uses a statistical band model with the strong line version of the Curtis- 
Godson (Godson, 1955) approximation; the water vapor dimer effect in the 8-13 
window region; line-by-line precalculation of GO 2 transmittance including funda- 
mental bands, hot bands, and isotopes; tables of ozone transmittance calculated 
by Dr. N. Scott (private communication); a special treatment of the nearby 
layer quadrature for increased accuracy in the numerical integration; and the 
effects of clouds. In order to calculate the incoming long wave flux at the 
10 mb level, the top of the model, zonally averaged climatological temperature'? 
at 1 mb and 5 mb are prescribed. In order to save computing time, the time 
step for the long-wave radiation is five hours for the simulations presented in 
this paper. 

Clouds are assumed to occur if and only if the model predicts cumulus con- 
vection (restricted to the lowest six layers) or large-scale saturation (in any 
layer). In keeping with the earlier discussion, no parameterization of subgrid- 
scale fractional cloudiness is attempted. In the control run described in this 
paper, both supersaturation and convective clouds are assumed to completely 
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fill a gcia box. In the experiment, superaaturation clouds are again assumed 
to completely fill a grid box, but convective clouds are ignored. For the 
solar radiation parameterization, the optical properties of the clouds are 
prescribed as shown in Table 1. For the terrestrial radiation parameterization, 
all clouds are assumed to be black bodies. 

In both runs, the modal was initialized with the observed conditions for 
15 June, 1979, and integrated for 45 days. The July mean results are obtained 
by averaging the output data through the last 31 days. Although history records 
were written only every 12 simulated hours, the radiation, cloudiness, and 
precipitation results presented in the next section ware accumulated between 
writes to the tape, so that there are no sampling errors in these results. 
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3. ROBUltS 

The most significant differences between the two simulations are in the 
radiation budget, as exiseeted, and in the temperature field. The feedback 
between cumulus convection and radiation also appears to be significant. We 
first show the cloud fields and the radiation budget; the effects on the other 
fields then become clear. 

Cloudines s. Tlie cloud field has been analysed to determine the global 
distribution and the aonal average of the July-mean frequency of cloudiness. 

The monthly mean cloud frequency is defined as the fraction of the time for 
which cloud occurs within the month. It has been obtained separately for cumulus 
clouds and supersaturation clouds. 

Pigs, la and 1b show the geographical distribution of the cumulus cloud 
frequency in layer 6 (about 500 mb), for control and experiment, respectively. 
Cumulus clouds occur somewhat less frequently over land in the experiment than 
in the control. This is true mainly for the deep clouds. The zonal-mean 
cumulus cloud frequency (Pig. 2) is reduced in the experiment everywhere over 
land by a few percent, except near BO^N and 30"S. Detailed examination of the 
model results shows that if we consider only deep cumulus clouds, the zonal-mean 
cloud frequency is reduced in the experiment not only over land, but also over 
ocean. 

When the model treats cumulus clouds as having 100% fractional cloud cover, 
there typically is strong long-wave cooling near the cloud top, and solar 
heating lower in the cloud. Radiative heating and cooling thus tend to steepen 
the lapse rate in the cloud layer, favoring stronger moist convection. This is 
why neglecting the radiative effects of cumulus cloudiness leads to reduced 
cumulus activity in the model. 




Helfand (1981) notod that the GLhS model ptoduffos *vjry intense cumulus 
oonvoction over land in summer. He showed that the simulated level of convective 
activity can be reduced by requiring that the relative humidity at the cloud 
base level exceed 95% as a precondition for convection. Our results show that 
the simulated convective activity over land in summer can also be reduced by 
neglecting the radiative effects of the cumulus clouds. Undoubtedly there are 
still more ways to reduce the simulated level of convective activity. 

For both the experiment and the control/ the supersaturation cloud frequency 
over the oceans is excessive in comparison with observations. Pig. 3 shows that 
much of the excessive supersaturation cloudiness occurs in the lowest model 
layer. The zonal-mean sup^y^aturation cloudiness (Pig. 4) is increased in the 
experiment over the ocean (but not over land). The increased supersaturation 
cloudiness tends to make up for the elimination of the cumulus cloudiness over 
the ocean. At the land points, there is a net reduction in cloudiness. 

^•2 Radiation . At the top of the atmosphere, the zonal-mean net radiation 
flux into the atmosphere (Fig. 5) is larger in the experiment, by about 20 W m~2, 
between 10-6Q®N. In the experiment, the local minimum of the radiation flux at 
about 26®N is due to increased oceanic supersaturation cloudiness. The results 
of the control run are in better agreement with observations obtained from 
satellite data of June 1974-February 1978 (Winston, et al . , 1979)^ However, a 
separate examination of the solar and thermal radiation streams indicates three 
characteristics of the radiation balance (Figs. 6 and 7). First, for each 
component of the radiation flux, the experiment shows better results than the 
control. The control underestimates both the outgoing thermal radiation and 
the incoming solar radiation more severely than the experiment, but these 
errors tend to cancel; as a result, the control produces a better simulation of 
the net radiation balance of the earth-atmosphere system. Second, the larger 




not incoming jra4'ltition fXux in tho ouinmtr hamispharc in bha expasimant/ raiativa 
bo bhab in bha conbrol, is chiefly duo bo a redvcbion in the albedo; the cloud 
albedo effect dominates the cloud greenhouse effect* Third/ both runs systema- 
tically underestimate the nob outgoing thermal radiation by about 40 W at 
all latitudes. We will come back to this point late»*. 

At the earth's surface/ the aonal-mean net radiation flux (Fig. 8) in the 
experiment is larger than that in the control/ especially near the ITCZ (approx- 
imately 10"N) and in the mid- latitude cyclone zone (30-60“^)* This is mainly 
due to reduced screening of solar radiation in the experiment. The difference 
reaches 15 W m"*^/ which is about 8% of the total flux. This difference occurs 
mainly over land/ where the model produces more cumulus clouds. In the experi- 
ment, the zonally averaged surface radiation flux over land increases, relative 
to the control run, by about 20% of the total flux. For reference, the observed 
surface net radiation flux is also shown (Schutz and Gates, 1972), although it 
is important to keep in mind the uncertainties of the observations. The exper- 
iment produces more solar radiation into the earth’s surface than the control 
in the summer hemisphere (Fig. 9), except near 26«N where the experiment has 
increased supersaturation cloudiness. Fig. 10 shows that the two simulations 
have about the same surface thermal radiation flux. 

Fig. 11 shows that the net surface radiation flux has increased in the 
experiment by about 20-60 W m"^ over Asia, North America, and central Africa. 
These locations of large differences coincide well with those of maximum fre- 
quency of cumulus clouds (Fig. 1). 

Herman (1980) suggested two possible explanations for the apparent failure 
of the GLAS model to emit enough thermal radiation to space. First, he noted 
that in the model cumulus clouds are unrealistically assumed to cover the entire 
grid area. In our experiment, we have completely removed the radiative effects 
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oC cumulus clouas, and in moat latitudinal belts/ the treaulfcing increase in the 
outgoing thermal radiation ia about 10 W (Pig. 7). h further inerenne of 
30 W m“2 is needed. 

Second, Herman pointed out that the model does not properly take into 
account the semi-branoparent x^^operties of thin, high cirrus clouds; these 
clouds are unrealistically assumed to bo black bodies, ox)aciue \.u the terrestrial 
radiation upwelllng from the earth's surface. In a separate eKx>eriment, shukla 
(personal communication, 1981) made all clouds in the top three layers of the 
model (above about 300 mb) transparent to both solar and terrestrial radiation, 
without changing the radiative effects of the cumulus clouds, which in any case 
do not penetrate the top throe layers. In Shukla* s exx»eriment, the net upward 
longwave radiation at the top of the atmosphere increased by about 30 W m“2 on 
aa average between SO^S and 30*N, and increased by about 5 W m’*2 ©n the average 
at middle latitudes. The tropical fluxes wore thus brought into reasonable 
agreement with observations, but the middle-latitude fluxes remained much too 
small. In contrast, the results of the present experiment (Fig. 7) show rela- 
tively little change in the tropics, but a noticeable Improvement at middle 
latitudes. 

Wo have found that part of the deficiency in the outgoing thermal flux to 
space is duo to the model's overprediction of low-level supersaturation clouds 
In both simulations, the July cloud frequency for supers^turation clouds in the 
lowest model layer over the ocean is twice that observed (see again Pig. 4). 

The simulated low stratiform cloudiness is about 20% more than observed, south 
of 30®N. These low-level clouds strongly omit thermal radiation downward to 
the earth's surface. For black-body emission, and assuming a 10“C temperature 
difference between the surface and the top of the low clouds, the excessive low- 
level clouds lead us to underestimate the simulated outgoing thermal radiation 
by about 12 W m-2. 
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3*3 Temporature and winds . Bocausa of the larger radiation input into the 
earth's ourface# the ground temperatures over Asia/ North America, and Africa 
are 3-6«C warmer in the experiment (Fig. 12). 0?his is in agreement with the 
results of Meloshko and Wetherald's (1981) study of the effect of a prescribed 
geographical distribution of clouds on climate. In their experiment, a reduc- 
tion in continental cloudiness results in warmer continents. As in their 
study, lower sea level pressures (Pig. 13) are obtained in our experiment over 
Asia, North America, and Africa. Over Australia, very few cumulus clouds 
occur; the warming in northern Australia is duo to a change of the large-scale 
low-level wind direction from southeasterly in the control to easterly in the 
experiment, which can be seen from the sea level pressure field. 

In the experiment, the atmosphere is generally warmer in the Northern 
(summer) Hemisphere (Fig. 14). This warming is due to the larger net radiation 
input at the top of the atmosphere (see again Pig. 5), which results from the 
fact that less solar radiation is reflected back to space in the experiment. 

The largest increase occurs in the upper troposphere^, near the tops of the 
deep cumulus clouds. Detailed examination of the output from the control simu- 
lation shows that the magnitude of the long-wave cooling near the top of the 
clouds is typically larger than the solar heating inside the clouds. Therefore, 
long-wave cooling dominates the solar heating, especially in the upper tropo- 
sphere. The experiment includes no radiative heating or cooling inside the 
cumulus clouds, and hence has warmer air in the upper troposphere. 

The higher temperature in the Northern Hemisphere in the experiment reduces 
the zonal-mean meridional temperature gradient, and therefore decreases the 
vertical shear of the zonal-mean zonal wind in the subtropics. Consequently, 

"I At the poles, there is only one grid point. Therefore, large differences 
between two simulations at the poles are not meaningful. 


’X2T 
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the zonal-mean zonal wind in the experiment is easterly throughout the whole 
subtropical troposphere (Fig. 15). 

3.4 Precipitation . Both simulations give about the same zonal-mean total 
precipitation (Fig. 16), except in the tropics and in the summer monsoon region 
near 26«N. The experiment produces increased convective precipitation in the 
coastal regions of the warm continents, near 26‘>N, as a result of the intensi- 
fication of monsoonal cirsulations associated with the warmer land in the experi 
ment. Meleshko and Wetherald (1981) reported a similar finding. This suggests 
that monsoonal cloudiness may serve, as a regulator of the monsoon circulations. 

A weak monsoon may produce relatively little cloudiness, allowing the continent 
to warm up, and thereby leading to an intensification of the monsoon. The 
extensive cloudiness associated with a strong monsoon may reduce the land-sea 
thermal contrast, thus tending to weaken the monsoon circulation. 

Fig. 17 shows the zonal-mean precipitation due to convective clouds only? 
here the control has a well-defined ITCZ but the experiment does not. By 
analyzing the convective precipitation over land and ocean, separately, we find 
that the control produces the tropical convection maximum at 10“N over both 
land and ocean, while the experiment produces the maximum over land at 10°N 
and a maximum over the ocean near the equator. As a result, the experiment 
does not show a well-defined tropical peak of the zonal-mean convective precipi- 


tation. 



4> Summary 

In the real atmosphere, convective clouds typically have much smaller 
horizontal extent, and hence much leas effect on the radiation balance, than v 

supersaturation clouds. For both types of clouds, the GLAS climate model i 

assigns 100% fractional cloud cover for a grid area to compute the radiation 
field. There is currently no method to parameterize the fractional cloudiness 
associated with deep convection. It may be better to assume zero cloudiness 
for convective clouds and 100% for supersaturation clouds. In this paper, we 
report a comparison between an experiment in which cumulus clouds are assumed 
to be radiatively inactive, and a control in which cumulus clouds are assumed 
to produce 100% fractional cloud cover. 

The control produces a closer agreement with the satellite-observed net 
radiation flux at the top of the atmosphere than the experiment. By separately 
analyzing the solar and thermal radiation fluxes at the top of the atmosphere, 
we find that the control underestimates both observed solar input and thermal 
output more than the experiment does, but by a cancellation of the larger 
errors, produces a better net radiation budget of the earth-atmosphere system. 

Over the whole globe, the zonal-mean thermal radiation flux from the earth's 
surface to the atmosphere and from the atmosphere to space are both systemati- 
cally underpredicted by about 40 W m"^ in both simulations, though the deficit 
is smaller in the experiment. These discrepancies partly result from an over- 
estimation of low-level supersaturation cilouds in the model. To significantly 
improve the model's simulation of cloudiness, it is therefore necessary not only 
to improve the parameterization of convective fractional cloudiness, but also 
to improve the parameterization of the turbulent transfer processes in the 
planetary boundary layer. An improved planetary boundary layer parameterization 
may prevent the overproduction of supersaturation clouds in the lowest layer. 


The comparison between two simulations shows a significant impact of the 
cumulus cloud -radiation interaction on model climatology. The cloud albedo 
effect dominates the greenhouse effect, and so neglecting the cumulus cloudiness 
produces more net radiation input into the earth -atmosphere system. As a con- 
sequence, the zonal-mean temperature is higher in the experiment. The warmer 
air occurs mainly over land, since the sea surface temperature is fixed in the 
model. It is clear that the climatic effect of the cumulus cloud-radiation 
interaction will be even larger in a coupled ocean -atmosphere climate model. 

In this light, efforts to parameterize subgrid-scale cumulus (and stratus) 
cloudiness are seen to be of great importance. 
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Table 1 . Prescribed cloud optical properties used in the GLAS solar radiation parameterization 
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Figure 1 j 

Figure 2 j 
F igure 3t 

Figure 4 j 

Figure 5s 

Figure 6s 
Figure 7; 

Figure 8i 

Figure 9s 
Figure 10s 

Figure 1 1 s 

Figure 12s 

Figure 13$ 
Figure 14s 
Figure 15s 
Figure 16s 
Figure 17s 


SimuXatod July convective cloud frequency for model layer 6 (about 
500 mb) for (a) the control and (b) the experiment. The contour 
interval is 20% and values greater than 40% are hatched. 

Simulated July sonal-menn convective cloud frequency averaged (a) 
over ocean and (b) over land. 

July zonal-mean supersaturation cloud frequency in the lowest model 
layer. The observations are from I^ondon (1957) and Van Loon et 
al. (1972), as shown in Meloshko and Wetherald (1981). 

Simulated July zonal-mean supersaturation cloud frequency averaged 
(a) over ocean and (b) over land. 

Zonally -averaged net total radiation flux fr<Mn space into the atmo- 
sphere. The observed data are from Winston et al. (1979). 

Same as Fig. 5, except for net solar radiation only. 

Same as Fig. 5, except for net terrestrial radiation only (positive 

upward to space). 

Zonally-averaged net total radiation flux into the earth's surface. 

The observed data are from Budyko (1963), as given by Schutz and 
Gates (1972). 

Same as Fig. 8, except for net solar radiation only. 

Same as Fig. 8, except for net terrestrial radiation only (positive 

upward from the earth’s surface). 

Global distribution of the difference (the experiment minus \*,he con- 
trol) in the net radiation flux into the earth's surface. The contour 
interval is 20 W m”^» Values greater than 20 W m“2 are hatched. 

Global distribution of the difference (the experiment minus the 
control) of the ground temperature. The contour interval is 2"C. 
Values greater than 2^C are hatched. 

Sea level pressure (a) of the control and (b) of the experiment. 

The contour interval is 4 mb. 

Difference (the experiment minus the control) of the zonal-mean 
temperature. 

Latitude-height section of the zonal-mean of the zonal wind compo- 
nent (a) for the control, and (b) for the experiment. 

Zonal-mean of the observed and simulated July mean total precipita- 
tion. The observations are from Jaegger (1976). 


Zonal mean of the simulated July mean convective precipitation. 
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